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Cuprates show multiphase and multiscale complexity that has hindered physicists search for the 
mechanism of high Tc for many years. Recently the interest has been addressed to a possible 
optimum inhomogeneity of dopants, defects and interstitials, and the structural scale invariance 
of dopants detected by scanning micro-x-ray diffraction has been reported to promote the critical 
temperature. In order to shed light on critical phenomena on granular materials, here we propose a 
stylized model capturing the essential characteristics of the superconducting-insulator transition of 
a highly dynamical, heterogeneous granular material: the Random Transverse Ising Model (RTIM) 
on Annealed Complex Network. We show that when the networks encode for high heterogeneity of 
the expected degrees described by a power law distribution, the critical temperature for the onset of 
the superconducting phase diverges to infinity as the power-law exponent 7 of the expected degree 
distribution is less than 3, i.e. 7 < 3. Moreover we investigate the case in which the critical state of 
the electronic background is triggered by an external parameter g that determines an exponential 
cutoff in the power law expected degree distribution characterized by an exponent 7. We find that 
for g = Qc the critical temperature for the superconducting-insulator transition has a maximum if 
7 > 3 and diverges if 7 < 3. 

PACS numbers: 89.75.-k, 89. 75. He, 05.30.Rt 



Large attention has been recently addressed to the 
effects that different topological properties may induce 
on the behavior of equilibrium and non-equilibrium pro- 
cesses defined on networks [l|, 0] . Heterogeneous degree 
distributions, small world and spectral properties, in par- 
ticular, have been recognized as responsible of novel types 
of phase transitions and universality classes in classical 
processes . For instance, scale- free networks present 
a complex critical behavior of the Ising model, of percola- 
tion and of spreading processes, that explicitly depends 
on the exponent of the power-law in the degree distri- 
butions [l|-|3[. On the other hand, the existence of non 
trivial spectral properties is crucial for the stability of 
synchronization processes and 0{n) models Never- 
theless there are no available results for quantum critical 
phenomena @ defined on complex networks. 

The aim of this paper is to investigate the role topo- 
logical aspect of the network on the Random Transverse 
Ising Model . In the last years there has been increasing 
interest in the understanding of the Random Transverse 
Ising Model This model that includes the dis- 

order of onsite energies, has been used to characterize 
the insulator-superconducting phase transition in gran- 
ular disordered superconductors 0] . In each grain of 
granular materials the superconducting order parameter 
is well defined and the grains are coupled to each other by 
the pair transfer term, so the physics is similar to the su- 
perconductivity in Josephson junction arrays The 



new results on a quenched Cayley tree have found the 
complex phase diagram of the model by the use the quan- 
tum cavity method The solution of the model 
is able to reproduce correctly the most important ex- 
perimental features on disordered granular films: direct 
superconductor-insulator transition^ctivated behavior 
close to the quantum critical point |12j in the insulating 



phase, strong dependence of the activation energy near 
the quantum critical point and the order parameter of the 
superconducting phase variations from site to site [13 1. 

The cuprate superconductors that keep the record of 
the highest critical temperature (160 K) of the phase 



diagram for all cuprate perovskite families jl4l | show a 
myriad of mysterious 'quantum matter' phenomena mak- 
ing these systems among the most intriguing puzzles in 
modern physics. In these doped Mott insulators the elec- 
trons form a poorly understood, highly collective quan- 
tum states. The strange electron matter of the cuprates 
strongly diverges from the standard weakly interacting 
quantum gas of conventional metals and superconduc- 
tors. Here the electron matter is in a 'quantum-critical 
state', where the electrons form collective patterns that 
look the same regardless of scale. This scale free topol- 
ogy is present both in space and in time scales, because 
the electrons are in a state of perpetual quantum motion 
IE EE EE]- ^ key particular feature of the anomalous 
electronic matter of cuprates is the heterogeneity of the 
interactions in different spots of the k space, going from 
weakly interacting electrons the so called " nodal points" 
to strongly interacting electrons at the " antinodal points" 
[l7| . These materials show multiphase complexity (struc- 
tural, magnetic and electronic) that has hindered physi- 
cists search for the mechanism of high Tc for many years 
since the observed phenomenology depends on the time 
and spatial sensitivity of the experimental probes. The 
heterogeneity of cuprate superconductors has been pro- 
posed to be an essential feature of high Tc mechanism 
jisl [lij and it has been proposed that above Tc the elec- 
tronic structure in the normal state is made of incoher- 
ent superconducting grains [I8, 2^, 21 1. In fact cuprates 
show a similar phenomenology for the essential features 
of superconductivity in disordered systems (i) the spot 



2 



to spot spatial variation on of the superconducting gap 
in cuprates detected by scanning tunneling microscope 
(STM) recording the static picture of states pinned at the 
dopant sites at the sample surface , and (ii) the recent 
characterization of the superconductor-insulator transi- 
tion by increasing continuously the charge density by gate 
voltage when the normal state resistance decreases below 
the value of the quantum resistance for pairs Rq = 6.5fcf2 

The focus is recently addressed on controlling oxygen 
defects both interstitials and vacancies to provide many 
avenues to control superconducting functionalities [24 1. 
In fact it has been shown recently that oxygen defects in 
heterostructures control other electronic states, including 
magnetism and ferroelectricity (25l |. a two-dimensional 
electron gas with universal sub-bands can arise from the 
oxygen vacancies confined in a surface layer at atomic 
limit in a cleaved SrTiOS crystals [26l - [28| and it could 
show superconductivity [2^. It has bee proposed that 
the defects clustering in cuprates might provide a form 
of optimal inhomogeneity for superconductivity [s^l • A 
fingerprint of structural scale invariance of dopants has 
been reported recently by Fratini et aL[3l| by detecting 
the structural scale invariance of dopants using scanning 
micro x-ray diffraction, a mixed fc-space and real-space 
probe, and confirmed by Poccia et al. [3^ using time 
resolved x-ray diffraction study of their nucleation and 
growth. The lattice complexity has been related with 
quantum criticality |33| . Both the magnitude and spatial 
distributions of the dopants distribution show power-law 
behavior with an exponential cutoff depending on the 
material and the doping, the unique fingerprint of scale 
invariance. Moreover the critical temperature Tc for su- 
perconductivity in these materials increases by a factor 
of 2.5 from 16K to 40K with the value of the exponen- 
tial cutoff changing by a factor of 3.8 [sH suggesting that 
the fractal background favors the superconductive phase. 
The scientific interest on this new scenario is rapidly 
growing since the control of complex structural organiza- 
tion of dopants and defects in new functional materials 
could allow us to manipulate granular superconductors 
with many venues for both science and technology. 

In this work we consider a possible scenario 
where at temperatures higher than Tc there are super- 
conducting grains without phase coherence (that we will 
call the insulator phase) and a superconducting coher- 
ent phase state below Tc (that we will call the supercon- 
ducting phase). We propose that the RTIM in annealed 
heterogeneous networks sheds light on the increase of the 
high critical temperature for an optimal heterogeneity in 
this scenario. Our model will describe the dynamical na- 
ture of networks that rapidly evolve in time by a rewiring 
of the links. Moreover a key element of our analysis is to 
study the role of heterogeneous degree distribution, fo- 
cusing in particular on scale- free degree distribution with 
an exponential cutoff which indicates the distance from 



a critical point which triggers the complexity of the elec- 
tron background. This will open a theoretical roadmap 
to mimic dynamics occurring on a complex fractal dis- 
order (31I I and will help to contribute in the general un- 
derstanding of how dynamical processes are affected by 
complex network topologies. We observe a rich interplay 
between network structure and quantum dynamical be- 
havior in annealed complex networks. In particular we 
have found that the critical temperature for supercon- 
ductivity is modulated by the topology of the underlying 
networks similarly to what happens in the classical Ising 
model on scale-free networks where it has been shown 
[5 that the phase diagram depends on the power-law ex- 
ponent of the degree distribution. Moreover the critical 
temperature depends on the power-law exponent on the 
degree distribution and on the exponential cutoff in the 
degree distribution that mimics the correlation between 
dynamical granular patches in cuprates. These results 
provide a new perspective showing how complexity can 
increase the critical temperature in an unconventional 
superconductor 24. 3^. 

Annealed Complex Networks - We consider networks 
of N nodes i = 1,...,A^. We assign to each node a 
hidden variable 6i from a p{9) distribution indicating the 
expected number of neighbors of a node. The probability 
(i,j) that two nodes are linked pij is given by 



Pij 



{e)N 



(1) 



In this ensemble the degree ki of a node i is a Poisson 
random variable with expected degree ki — Oi . Therefore 
we will have 



= (fc(fc-l)). 



(2) 



Random Transverse Ising model on annealed complex 
networks - We consider a system of spin variables , for 
i = 1, . . . , Af, defined on the nodes of a given annealed 
network with link probability given by the matrix p and 
adjacency matrix a. The Random Traverse Ising Model 
is defined as in [1, Q as 



(3) 



This Hamiltonian is a simplification respect to the XY 
model Hamiltonian proposed by Ma and Lee 01 to de- 
scribe the superconducting-insulator phase transition but 
to the leading order the equation for the order parameter 
is the same, as widely discussed in d, The Hamilto- 
nian describes the superconducting-insulator phase tran- 
sition as a ferromagnetic spin 1/2 spin system in a trans- 
verse field. We propose to use this Hamiltonian to de- 
scribe in a granular superconductor the transition from 
a phase of superconducting grains with no phase coher- 
ence (called insulator for granular superconductors) that 
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we propose to correspond with the electronic matter in 
cuprates above to the low temperature superconduct- 
ing phase with phase coherence. 

The spins Oi in Eq. ([3]) indicate occupied or unoc- 
cupied states by a Cooper pair or a localized pair; the 
parameter J indicates the couplings between neighbor- 
ing spins, ei are quenched values of on-site energy and h 
is an external auxiliary field. To mimic the randomness 
of on-site energy we draw the variables from a p{e) 
distribution with a finite support. Finally in this model 
the superconducting phase corresponds to the existence 
of a spontaneous magnetization in the z direction. The 
partition function for this problem is given by 

Z = Tr e-f^" (4) 

with the Hamiltonian given by Eq. (j3]) where in order to 
account for the dynamics nature of the annealed graph 
we have substituted the adjacency matrix aij in H with 
the matrix pij given by Eq. [T] In order to evaluate the 
partition function we apply the Suzuki- Trotter decompo- 
sition [s^l in a number Ns of Suzuki- Trotter slices. This 
formula expresses the exponential of the sum between two 
operators in terms of a limit of a product of exponentials, 
i.e. 



„A+B 



= lim [e 



A/N^^B/N,-\N, 



(5) 



Therefore we have in the limit Ns oo that the partition 
function Z can be written as 

Z = Tv e-l'" = Tr (e-^^/^=e-'3^'^'*'/^0'^° , (6) 



where E is given by 
J 



(7) 



where are the spin operators. In order to perform 
the calculation of the partition function Z we consider 
for each spin the sequence ct, = {a}, . . . af"} where each 
spin represents the spin i in the Suzuki- Trotter slice 
a. The partition function is then defined as 



N 



I3h v^iVs 
5 JV. 2^0 = 1 ' 



{2.i}i=l.....N ' = 1 



13. J 



where we have indicated with 



(8) 



(9) 



In order to disentangle the quadratic terms we use Ns 
Hubbard-Stratonovich transformations 



N,/2 



Z . r«V^-'" ^^exph^E(^") 



V 27rA^s / 

exp[Nj2pi&) I rfep(e)lnTrJ| 



2N:_ 



where VfS 



The free energy / 
— limAT^oo limAT^^oo In .Z can be evaluated at the sxa- 
tionary saddle point which is cyclically invariant. There- 
fore we get 

/ = mfs-^-S^ 

-^^^(61) I dep{e) In (2 cosh {p^J{h + JSOy + e 

P g J 

where the value of S which minimizes the free energy is 
given by the saddle point equation 



jse + h 



(10) 



Finally the magnetizations along the axis a; and z can be 
calculated by evaluating 

Tr afe-l^^ 



Z 

Tr cjfe-^" 



Z 



(11) 



Performing these calculations we get 



jse + h 



^{jse + hf 
^{jse + hY 



Therefore the magnetizations rrig ^ and rrvg ^ depend on 
the value 9 of the expected degree of the node and on the 
on-site energy e. 

These magnetization curves are plotted as a function 
of 9 and e in Figure [1] for some parameters values. The 
order parameter for the superconducting-insulator phase 
transition is S given by Eq. ([T0| . 

From the self-consistent equation determining the or- 
der parameter for the transition it is immediate to show 
that the superconducting-insulator phase transition oc- 
curs for ft, = at 



(0) 



J 



JcW) 



(12) 



which implies that for {9'^) ^ 00 then /? — !> for any fixed 
value of the coupling J > , and the critical temperature 
for the paramagnetic ferromagnetic phase transition Tc 
diverges. This implies that on annealed scale free net- 
works with p(9) oc 9^^^ and 7 < 3 the random transverse 
Ising model is always in the superconducting phase. 

In order to mimic the fractal background present in 
cuprates [3l[ we assume that the expected degree distri- 
bution of the network is given by 



p{9) 



(13) 



FIG. 1: (Color online) Magnetization rrig ,, (cyan top surface) 
rrig £ (green lower surface) as a function of and e for T — 10, 
expected degree distribution p{6) oc k^'' and 7 = 2.5, and 
distribution of on-site energies p(e) = 0.5 with e — (—1, 1). 

where ^ eventually becomes much larger than the lattice 
coordination number and diverges as a function of an 
external parameter g (mimicking doping of the cuprate 
or strain of the lattice) with a critical exponent v, i.e. we 
assume 

e = lff-ffc|-^ (14) 

In figure [2] we show the behavior of the critical tempera- 
ture as a function of the external parameter g below and 
above the critical point gc- It is shown that as 7 > 3 the 
critical temperature of the RTIM has a maximum at Tc 
while for 7 < 3 it diverges at g = g^ 

Conclusions - In this paper we have investigated 
how topological effects and heterogeneity in the de- 
gree distribution can affect the critical behavior of the 
superconductor-insulator transition in annealed complex 
networks. In particular, we have shown how topologi- 
cal effects change significantly the phase diagram of the 
critical phenomena. Our analysis points out that a finite 
second moment of the degree distribution determines the 
critical lines of the transition while when this second mo- 
ment diverges the critical temperature for the ferromag- 
netic phase diverges. In the future we plan to study how 
much this scenario changes if we consider quenched com- 
plex networks by applying the quantum cavity method. 
Finally our results open new perspectives for the com- 
prehension of critical phenomena in complex materials 
when a structural phase transition driven by some exter- 
nal parameter modulates with an exponential cutoff 
the scale-free heterogeneities in the material. We found 
that in this case two scenarios are possible: either the su- 
perconducting temperature reaches a maximum at g = gc 
or the superconducting temperature diverges at g ^ gc 
depending on the value of the power-law exponent 7 of 
the degree distribution of the annealed network. We 
hope that this work will be of help for the design of new 
complex granular high temperature superconductors that 
could eventually overcome the limit of 160K in cuprate 
superconductors. 



FIG. 2: (Color online) Left panel: The critical temperature Tc 
for the Random Transverse Ising Model with J — 1 p(e) = 0.5 
with e £ ( — 1, 1) as a functional of the external parameter g 
tuning the exponential cutoff of the expected degree distribu- 
tion p{e) (X e-^e-"^^ with C = Iff - gc\''",9c = 1 and = 1. 
The critical temperature Tc for the superconductor insulator 
transition has a maximum for g — gc — lif'y>3 and diverges 
at (? = (?c is 7 < 3. Right panel: critical temperature Tc for 
the Random Tranverse Ising Model with J = 1 p(e) = 0.5 
with e £ ( — 1, 1) as a functional of the exponential cutoff in 
the expected degree distribution for different values of the 
power-law exponent 7. 
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